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Abstract 

A consistent implementation of quantum gravity is expected to change the fa- 
mihar notions of space, time and the propagation of matter in drastic ways. This 
wiU have consequences on very small scales, but also gives rise to correction terms 
in evolution equations of modes relevant for observations. In particular, the evolu- 
tion of inhomogeneities in the very early universe should be affected. In this paper 
consistent evolution equations for gauge-invariant perturbations in the presence of 
inverse triad corrections of loop quantum gravity are derived. Some immediate ef- 
fects are pointed out, for instance concerning conservation of power on large scales 
and non-adiabaticity. It is also emphasized that several critical corrections can only 
be seen to arise in a fully consistent treatment where the gauge freedom of canonical 
gravity is not fixed before implementing quantum corrections. In particular, metric 
modes must be allowed to be inhomogeneous: it is not consistent to assume only mat- 
ter inhomogeneities on a quantum corrected homogeneous background geometry. In 
this way, stringent consistency conditions arise for possible quantization ambiguities 
which will eventually be further constrained observationally. 



1 



1 Introduction 



General relativity describes the structure and dynamics of space-time by Einstein's equa- 
tion for the space-time metric, which gives rise to a wide range of phenomena in cos- 
mology and astrophysics. In several regimes, especially at high densities, one expects 
quantum gravity to be crucial and to provide important correction terms to the classi- 
cal equations. By now, canonical quantum gravity in the loop approach has progressed 
to the extent that perturbative calculations for the behavior of inhomogeneities around 
Friedmann-Robertson-Walker space-times can be performed at an effective level. This 
provides cosmological applications, but by addressing the anomaly problem it also sheds 
light on fundamental aspects such as quantum space-time structure. 

The classical equations of motion constitute an overdetermined set, whose consistency 
is ensured by general covariance. When attempting to include correction terms in these 
equations as they may be suggested by quantum gravity, the consistency conditions must 
be respected. Only an anomaly-free quantization, where consistency resulting from general 
covariance is preserved, can lead to quantum equations of motion with the correct number 
and behavior of degrees of freedom. 

This issue becomes pressing already at the level of linearized equations as they are used 
in cosmology of the early universe. Once inhomogeneities are included as perturbations 
around an expanding Friedmann-Robertson-Walker space-time, there are more equations 
than unknowns which requires consistent forms of all terms in the equations. Without 
inhomogeneities, consistency is automatically satisfied: there is a single constraint which 
is always preserved by the evolution it generates. Quantum corrections to homogeneous 
models can therefore be implemented easily as it has been done often to suggest diverse 
effects and scenarios p. However, just putting quantum corrected solutions for a homo- 
geneous background into classical perturbation equations in general results in inconsistent 
equations: the corrected background equations can no longer be compatible with all the 
terms in the perturbation equations. Thus, a direct treatment of inhomogeneities and the 
corrections they acquire in a quantum theory of gravity is required. 

Perturbations around homogeneous models then test their robustness, demonstrate 
whether a particular form of quantum corrections can be realized in a generally covariant 
way, and provide consistent sets of equations of motion whose solutions can be analyzed for 
the phenomenology and potentially observable effects they imply. From the perspective of 
canonical quantum gravity, the consistency issue of effective equations has been described 
in [2], and it has been demonstrated that there is a correction expected from loop quantum 
gravity [31 lU E] which non-trivially changes the classical equations in a consistent way at 
a perturbative level. (This correction results from inverse triad operators in Hamiltonians 
[6l[7j. Other expected corrections from loop quantum gravity, such as holonomy corrections, 
have not yet resulted in consistent equations outside the reduced setting of homogeneous 
models.) In this paper, we derive the corresponding gauge invariant perturbations and the 
equations of motion they satisfy. 

Earlier work [H |9] had already led to quantum corrections to Einstein's equation 
governing linear cosmological perturbations. As a result, enhancement effects of quan- 
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turn corrections during long cosmic evolution times were suggested based on the observa- 
tion that super-horizon curvature perturbations were not preserved, unlike classically, but 
had a growing mode. Other forms of perturbation equations were used, for instance, in 
[TOl [TTl IT2I [131 lEy. Such effects may lead to observable imprints in the cosmic background 
radiation even though individual quantum-gravitational correction terms at sub-Planckian 
densities are small. The consideration, however, was restricted to the scalar mode in the 
longitudinal gauge, i.e. diagonal metric perturbations. While classically this procedure 
proves to be equivalent to a non-gauge fixed derivation, it leads to inconsistencies at the 
effective level in the presence of quantum corrections. Specifically, the effective equations 
resulting from the quantization of a gauge-fixed system in general are incompatible with 
each other, and gauge-fixing eliminates the freedom required to see systematically how the 
terms of consistent equations must be arranged. A non-gauge fixed treatment is thus neces- 
sary to evaluate all consistency conditions and to determine the gauge-invariant equation of 
motion for the curvature perturbation. (Similarly, if one uses only matter perturbations on 
a homogeneous gravitational background one is implicitly using a gauge-fixed treatment. 
This is in general inconsistent if no care is taken concerning the specific correction terms 
and the meaning of matter perturbations in relation to gauge-invariant quantities. Some- 
times a "separate universe" picture [151 [13 [U] is used, arguing that at least for large scale 
modes quantum corrections to the homogeneous background equations should be sufficient 
to determine the evolution of inhomogeneities. But also here, as we will see in examples 
later, not all features will be visible based solely on homogeneous models without a full 
consistency analysis once quantum gravity corrections are included.) 

A simple counting shows that the three independent functions describing scalar (gravity 
and matter) perturbations are subject to five equations of motion (see the Table [1]). In 
general, such a system would be over-determined unless the equations are not independent. 
In gravity, they can be split into two types: i) evolution equations of second order in time 
derivatives and ii) constraint equations of lower order. Constraints restrict the initial 
data, and if they are preserved under evolution the system of equations is consistent. 
This is guaranteed automatically when equations of motion are obtained variationally 
from a covariant action, and is therefore satisfied for the classical equations. If quantum 
corrections are derived in a Hamiltonian approach, however, consistency is ensured only if 
the quantization is anomaly-free. The consistency of the resulting equations is thus tightly 
related to the closure of the constraint algebra. While the algebra of the constraints 
obtained from the classical Einstein-Hilbert action by a Legendre transformation is closed 
or, in Dirac's terms the Hamiltonian and diffeomorphism constraints are of first class, 
this property may not sustain quantization. A consistent gauge invariant formulation of 
quantum-corrected equations of motion is possible only if the quantization is anomaly free, 
i.e. if the constraints remain first class. 

As shown in p!9| [20]. standard loop quantization under very mild assumptions leads to 
a non-anomalous constraint algebra for vector and tensor modes. (At the linear level mode 
decomposition does not interfere with quantization, and quantum corrections to scalar, 
vector and tensor modes can be studied independently.) In [2] it was analyzed what types 
of (non-anomalous) quantum corrections are allowed for the scalar mode, obtaining the 
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anomaly-freedom conditions ensuring a first class system. Once closure of the constraint 
algebra is provided, the formulation of gauge invariant equations of motion becomes pos- 
sible, as developed in this paper. 

We start with reviewing the correspondence between the canonical and covariant equa- 
tions of motion, then derive the gauge invariant variables and finally obtain the gauge 
invariant quantum corrected linear Einstein equations. As we will see, consistency requires 
certain features of the corrected perturbation equations and of the gauge invariant vari- 
ables, which could not be seen in gauge-fixed formulations. Several immediate consequences 
are discussed in the Sec. |5]and further in the conclusions, which also exhibit the final quan- 
tum corrected perturbation equations. Readers interested primarily in applications may 
turn directly to these sections. 

2 Basic variables and equations 

When a classical theory is quantized, the choice of basic variables often matters. While 
there are many equivalent formulations of classical physics, all related to each other by 
canonical transformations, such maps are rarely implementable as exact unitary transfor- 
mations when quantized. This gives rise to different inequivalent quantizations of the same 
classical theory, and it can even prevent one from constructing a quantization in a particu- 
lar classical formulation of the theory: there may be no Hilbert space representation where 
a certain choice of basic classical phase space variables will become well-defined operators. 

In loop quantum gravity [31 HI |5], the principle of background independence, which 
requires that well-defined operators do not refer to a metric other than the physical one 
to be turned into operators, distinguishes a special class of basic variables. In field the- 
ories such as general relativity, it is not field values at single points which can become 
well-defined operators, but only "smeared" versions obtained after integrating over spatial 
regions. Such integrations ensure that the operator- valued distributions, which field values 
would correspond to, become well-defined operators which can be multiplied to construct 
composite operators from them. Since the physical fields of canonical general relativity 
are spatial tensor fields, they cannot directly be integrated in a coordinate independent 
manner. Moreover, integration measures are not provided automatically because only the 
physical metric could be used on a curved manifold, but this metric itself is being turned 
into an operator. If one uses connection variables and densitized vector fields as canonical 
objects, however, their transformation properties ensure that they can be integrated over 
curves and surfaces, respectively, without requiring any additional integration measure. 
The resulting holonomies and fluxes then become well-defined operators in the quantum 
representation underlying loop quantum gravity [21j . 

This representation has characteristic properties which are implied by the choice of 
basic fields and their smearing. In particular, operators for spatial geometry such as the 
fluxes themselves or areas and volumes acquire discrete spectra [221 [221 121] • This, in turn, 
determines how these basic operators can appear in composite ones such as Hamiltonians 
[251 El Ej- For instance, flux operators having discrete spectra containing the eigenvalue 
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zero do not possess densely defined inverse operators. Since inverse triad components 
appear in Hamiltonians, the lack of a direct quantization entails quantum corrections in 
any effective Hamiltonian (see e.g. [26]), which will then also change the corresponding 
evolution as well as gauge properties. In this paper, we derive such equations precisely for 
corrections resulting from inverse triad components. 



2.1 Perturbed variables 

To do so, we perform the perturbation analysis of inhomogeneities in the basic variables 
underlying loop quantum gravity, such that we will be using primarily a densitized triad Ef 
instead of the spatial metric qat (satisfying E^E^ = q"-'^ det q) . Moreover, in this canonical 
setting the remaining components and A^" of the space-time metric 

ds^ = -N^dt^ + gafe(dx'^ + N^dt) {dx^ + N^dt) (1) 

will not be dynamical but play the role of Lagrange multipliers of constraints. In fact, 
their time derivatives do not appear in the Einstein-Hilbert action, which can be written 
in the canonical form 



5eH = jdt j d'xk^Et^ - Ggrav[Ai - I^grav[iVl " ^^grav [iV] 

where i^"* is conjugate to -Ef , 



(2) 



{Kl{x),E^{y)} = 87cG6l6}6\x,y) , (3) 

and related to extrinsic curvature Kab by K^^ = KabE^ / ^^/\'det^. The remaining terms 
are the diffeomorphism constraint 

and the Hamiltonian constraint 

1 f E°'E^ / \ 

- L '''''' 7fm " ^ ^^"'"^ ■ 

Here, 

F^,, = 2dUT + -fK)\ + 6./(r + ^K)\{T + ^K)i 

is the curvature of the Ashtekar-Barbero connection [271 [28] A\ = T^^ + 1^1^ defined in 
terms of the spin connection 



(6) 
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and 7 is the Barbero-Immirzi parameter [281 129] • These variables also appear in the Gauss 
constraint 

in (j2l). This constraint will be solved explicitly by our parameterization of variables at the 
linear level, and its gauge will be fixed by a background triad. We can thus ignore it from 
now on. 

For a perturbed metric of the form 

ds^ = a'ir]) (-(1 + 20)dr/2 + 29„5d?7da;" + ((1 - 2^)5,^ + 2dad^,E)dx''dJ') , (7) 

as it describes general scalar perturbations {(f),ilj, E, B) around spatially fiat Friedmann- 
Robertson-Walker models in a general gauge, the background and perturbed triad in 

Et = E^ + 6Et (8) 

are given by 

E^=p6t^a'6t , 6Et = -2p^6t+p{6tA-d''d.)E, (9) 

where A is the Laplace operator on fiat space, indices run from 1 to 3, a is the (background) 
scale factor and ip ^ind E describe the spatial part of the perturbed metric. (We use 
standard notations where E is one of the scalar modes distinct from the full densitized 
triad E^. The latter will always be written with indices such that no confusion should 
arise.) The other two scalar metric perturbations and B are related to the perturbed 
lapse function and shift vector respectively 

6N = N(j) , N" = d^B (10) 

and will enter the extrinsic curvature components, also perturbed as 

K = K + 6K = k6i + 6K. (11) 

This splitting, with the condition that SE"- and 6Kl^ do not have homogeneous modes in 
order to avoid double counting, results in Poisson brackets 

{kp} = ^ , {6K{x),6E';{y)} = 8nG6}6'j'{x-y). (12) 

The homogeneous mode is defined by 

where we integrate over a bounded region of coordinate size Vq = J d'^x which could be 
over the whole space if it is compact, or a sufficiently large region encompassing all the 
scales of perturbations of interest. Although Vq, which depends on coordinates as well as 
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the choice we make for the integration region, enters the definition of variables and their 
Poisson structure, particular correction terms for observables will not depend on its value. 
Using the homogeneous modes, perturbations are defined by (JHl) and (fTTl) . 

The specific form of and 6KI in relation to time derivatives of p and 6E^, analogously 
to the triad fields in ([9]), follows from the equations of motion. Before deriving these 
relations we thus introduce the quantum corrections we consider because they have a 
bearing on the form of components of /T*. 



2.2 Quantum corrections 

The Hamiltonian constraint ([5]) contains a factor of an inverse determinant of the densitized 
triad. This inverse cannot be quantized directly because the integrated determinant itself is 
quantized to an operator with zero in the discrete spectrum, precluding the existence of an 
inverse operator. Nevertheless, well-defined operators quantizing ([5]), including the inverse 
triad, exist [H [30] . However, the behavior of expectation values of the operators differs on 
small length scales from the classical behavior even in semiclassical states, which implies 
the presence of a correction function a multiplying the Hamiltonian density. This function 
must be scalar (of density weight zero) to ensure the proper behavior of the integral, and it 
can depend functionally on all the phase space variables in possibly non-local ways. Such 
a general dependence would make an analysis of the constraint algebra and of equations 
of motion intractable, and so we have organized the calculations in [2j by first assuming 
a primary correction function a(-Ef) which depends only on the triad and does so only in 
algebraic form. By itself, this does not produce anomaly-free quantizations, which however 
do exist if additional counter-terms are added containing new correction functions whose 
relation to the primary correction is fixed by anomaly-cancellation. These extra terms can 
be interpreted as arising from a more complicated dependence of a on all the phase space 
variables, which is derived systematically by this process. 

Specifically, the quantum corrected Hamiltonian constraint derived in [2] can conve- 
niently be written as 

= H^,.JN] + H^^J6N] + H^^JN] + h2^,,[6N], (14) 
where the gravitational part is expanded by powers of inhomogeneities SEf, SK^^ and 6N 



as 



//Q [AT] := [ d^xN + «(2)7^Q(0) + ^^Q(2)i ^ 

IGttG J 

HiJSN] ■■= d'x5N [an^W^ , (15) 



and the matter Hamiltonian reads 



H2..t..Am = j d'xdNlm^^^'^ + nf^] . (le) 



7 



In all expansions, a bar is used to denote background quantities while superscripts indicate 
the inhomogeneous order. The Hamiltonian densities in the expansions are given by [2] 

VP VP 
2k 

^ ^SKi6K',6l6^ - Vpmi6^r - -^6E^6Ki (17) 
for gravity and by 

Hg^'^ = (18) 



n2 



= ^(l + ^75)v^5'^U5<^5,5<^ 



^Q(2) ^ ^3/2 



1 , 6j.5E^ i^i^E'^? 6^5^JE^6E? 



for matter. 

The p-dependent functions a, v and a are primary correction functions whose origin is 
the presence of inverse triad operators in a constraint operator. Their form can be com- 
puted in isotropic models [SOj [311 [32] or with certain gauge assumptions for inhomogeneous 
states [331 [26]. Classically, we have a = z/ = a = 1, while there can be strong deviations 
from this value for small values of elementary flux variables which quantize the densitized 
triad. This deep quantum regime is difficult to control, however, and the derivations in 
[2] of an anomaly-free constraint algebra are valid for primary correction functions of the 
form 

/ «2 \ 

a(a) = l + cj^j +■■■ (19) 

which are perturbative in the Planck length £p = VGh, i.e. Ua > 0. Explicit values for 
coefficients c^, c^, and c^, which are generically positive such that a{a) > 1 in perturbative 
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regimes, as well as the exponents Ha, Uy and can be derived from specific quantizations, 
but they are subject to quantization ambiguities^ One purpose of deriving anomaly-free 
versions of the constraints is to provide consistency conditions among some of these values, 
fixing some quantization ambiguities. 

For an anomaly-free quantization in a gauge-independent manner, the presence of these 
primary correction functions requires counter-terms with coefficients /, g and h as well as 
fi and Qi which also depend on p in a way fixed by anomaly cancellation conditions. For 
the situation under consideration where the matter sector consists of a scalar field with a 
non-trivial potential, we have 

VP = -— (20) 
a 

9 = -2/ (21) 
/i = /-^ (22) 



and 

a' _ az/(2) 



Here and in what follows, primes denote derivatives by p. Moreover, we have 

oi^ = ua. (24) 

Other consistency conditions will be recalled later from [2] (also discussed in Appendix [B|) 
whenever they are being used. Classically, all counter-terms vanish, e.g. f = fi = g = 0. 
With the consistency conditions the system of corrected constraints is anomaly-free to the 
perturbative orders considered, which is linear in inhomogeneities (requiring second order 
expansions of the constraints which generate linear equations of motion) as well as leading 
order in the corrections of (fT9l) . The latter assumption of perturbativity implies that we 
ignore terms such as {a — 1)^, (pda/dp)"^ or compared to a — 1. 

With the corrected Hamiltonian, we can derive the equations of motion it generates. 
From p = {p, H^^^^[N]}, for instance, we obtain the background part K"* = kS^ of extrinsic 
curvature where k is related to the conformal Hubble parameter by 

ak = jr=^. (25) 
2p 



^Even the isotropic quantization used for the background evolution is subject to quantization ambigu- 
ities. Uniqueness results of the quantum dynamics can be obtained only based on ad-hoc assumptions, 
and they sometimes occur as a result of incorrect implementations of quantization schemes. The source 
of ambiguities is the representation of operators, such as inverse triad operators used here, but also the 
underlying refinement behavior of a discrete state underlying the quantum evolution |341 |35j : see also 
and the appendix of 2 . Both ingredients combine to determine the values of and ■ Note also 
that Cq, when correctly derived using lattice refinement, is coordinate dependent in such a way that the 
combination with the scale factor in pop is scaling independent. 
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The choice of the background lapse function N = a, used to derive (l25l) . corresponds 
to the conformal time t] whose derivative we denote by a dot. In general, the total time 
derivative of an arbitrary phase space function is given by its Poisson bracket with H'^[N] + 
DIN"] parameterized by the total lapse = N+6N and shift A^" = N^+SN". Nonetheless, 
for a background quantity the Poisson bracket above, using only H^[N], coincides with 
the conformal (background) time derivative up to the second perturbative order. 

Similarly, the form of the perturbation SK^^ can be deduced from Hamilton's equation 
for . Namely, using 

6Et = {6Ef, H^[N] + D[N'']} (26) 
along with and (^S^, we obtain 



cxSKl = -61 



^/^ + ^(V' + 0(1 + /)) +dad' J^E-{B-E) , (27) 



where the counter-term f{j>) appears. 

Matter is represented by a scalar field = (p + 5(p with potential V{(p) and its conjugate 
momentum tt = vr + dn. As before, we use the equations of motion 

= {ip,H'^[N]+D[N'']}, 5iP = {5ip,H'^[N]+D[N'']} (28) 

to express the field momentum as 

TT = 571 = ^ {{5^ - ^(1 + /i)0) (1 - 9i) + ^^f) ■ (29) 

Before proceeding to gauge transformations in the next section we note the relation 
between the canonical and covariant equations of motion, summarized in Table [H There 
are three background equations, only two of which are independent, for the two unknown 
functions: scale factor a{ri) and matter scalar field '^{rj) depending on the conformal time 
7]. Those are the Friedmann, Raychaudhuri and Klein-Gordon equations [26] 

= (I (30) 

^^^^(l.f)-4.«|,^(:-£|). (31) 

<^ + 2^^p - + upV,^{gy) = 0, (32) 

where the prime indicates a derivative with respect to p. These equations are listed in the 
left column and correspond to the background Hamiltonian constraint and two pairs of 
dynamical (Hamilton's) equations. Each pair of the first order equations is to be combined 
into a single second order equation. 

More generally, covariant equations which are less than second order (with respect to the 
conformal time derivative) correspond to constraint equations in the canonical formalism. 
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Covariant Equations 


Canonical Equations 


Background Friedmann 


Background Hamiltonian Constraint 


Background Raychaudhuri 


k &: p 


Background Klein- Gordon 




Perturbed Einstein ^ 


Perturbed Hamiltonian Constraint 


Perturbed Einstein ^ 


Perturbed Diffeomorphism Constraint 


Perturbed Einstein | 


6K k 6E 


Perturbed Klein-Gordon 


5ip & 5ti 



Table 1: Table of correspondence between the background and perturbed canonical equa- 
tions and the covariant equations of general relativity. In the covariant column, the sub- 
scripts 'S' and 'T' stand for spatial and temporal components respectively. In the canonical 
framework, equations are of two types: constraint equations and dynamical (Hamilton's) 
equations for the time derivatives of canonical pairs. Note that both 'SS'- and ^5K & 
'-equations are tensorial. Indices of 5Kl^ and 5E°; have been suppressed for simplicity. 



They should be viewed as restrictions on the initial conditions. As mentioned earlier, for 
consistent and unambiguous (gauge-invariant) evolution such constraints must be preserved 
by the dynamical equations. In the canonical language, this property of constraints is 
translated into the requirement of closure of the constraint algebra, as it is analyzed for 
the corrected constraints in [2J. As a result, relations such as (1201) . (1221) . (I2T]) and (!23|) 
between the correction functions have to be satisfied as conditions for higher order terms 
of primary correction functions. 

3 Gauge transformations 

In classical relativity, it is the Lie derivative which provides the form of gauge transfor- 
mations of the fields such as metric components, corresponding to changes of coordinates. 
After quantization it is no longer clear what the analog of these gauge transformations 
should be, given that the underlying space-time notion would have to be determined from 
the quantum theory itself. In most approaches to quantum gravity, one does not expect 
the fundamental space-time picture to be described by a smooth manifold. Here, an ad- 
vantage of the canonical formulation is that gauge transformations are directly generated 
as Poisson brackets of the fields with the constraints. Classically, this reproduces the for- 
mulas obtained by Lie derivatives, and it can directly be extended to canonical quantum 
gravity capturing changes to the quantum space-time structure. With corrections (fT4l) to 
the classical constraints ([5]), it is not only equations of motion but also the form of gauge 
transformations which changes. Thus, gauge invariant combinations of the perturbations 
take different forms than they do classically. 

However, as we saw not all the space-time metric components are dynamical phase 
space variables, and only the gauge transformations for the spatial metric, or E°; and i^'*. 



11 



will be determined straightforwardly. In this section we derive these transformations and 
show how also the transformations of the remaining components N and N"-, or and B 
in the scalar perturbations related to the Lagrange multipliers, can be obtained. 

3.1 Classical gauge-invariant variables 

It is instructive to introduce the canonical derivations and required notions in the classical 
case first, after which we will directly extend the expressions to those including quantum 
corrections. (Canonical treatments of classical perturbations have also been discussed in 
[371I3HII39].) 

In the covariant formulation, gauge transformations constitute local infinitesimal coor- 
dinate transformations 

xi' ^ = x" + 

generated by vector fields In a perturbative setting, the infinitesimal field ^'^ will be 
treated as a first order perturbation. Under this coordinate transformation any tensor field 
receives a correction equal to its Lie derivative along The part of the transformation 
relevant for the scalar mode can be parameterized by two scalar functions and ^ such 
that 

where a indicates a spatial direction 1, 2 or 3. 

There are four (spatial) scalar perturbations in a space-time metric, (p, ip, B and E 
as they appear in the line element ([7]) in conformal time 77. These perturbations are 
subject to two independent gauge transformations by C,^ and ^. We now briefly recall how 
these transformations follow from changes of coordinates to verify later that the classical 
canonical transformations produce the correct form. 

If only is non-zero, the coordinate transformation changes 77 to 77 -|- which for d?7^ 
implies, to first order in 

d(r/ + = dr]'^ + 2i^<lrf + 2^°,^ dr/dx" 

and a(?7)^ changes to a{riY{l + 2a^^ / a). Inserting this in ([7]), we read off the transformation 
formulas 

(/) (/) + ^° + , ip^ip-J^i"^ , B^B~i^ , E^E (33) 

where = a/a as in fl25l) . 

If only ^ is non-zero, dx^da;'* changes to 

d(a;„ + e,a )d(a;'^ + e," ) = dx.dx" + 2^ dr/dx" + 2Ub dx-^dx" 

which yields 

(/)^<^ , ip^^ , B^ B + i , E^ E + ^. (34) 
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We thus see that B—E is invariant under ^-transformations (or spatial diffeomorphisms) 
and changes to B — E — C,^ under ^^-transformations. Thus, the Bardeen variables [IQ] 



^ = ^ + J^{B - E) + (B - E)- and = ip - J^{B - E) (35) 

are gauge invariant. For a scalar field (p, the only change is under and given by 6(f + (f^^. 
Here, 

= _ E) (36) 

is gauge invariant. 

In the canonical formulation, gauge transformations are generated by the Hamiltonian 
and diffeomorphism constraints. The corresponding lapse function and shift vector to be 
inserted are also first order perturbations related to the infinitesimal vector field ^'^ via 

SN = N^o, SN" = (9"e, (37) 

which follow from the metric decomposition ([1]) in terms of the spatial metric Qab, lapse 
function and shift vector A^". With first order smearing functions, the gauge-generating 
constraints are at least of second perturbative order. From now on we will denote the 
gauge transformations of a phase space quantity X as 

^ {X,H('^mo]} + {X,D^'^[d'^^]}, (38) 

where 

H^^)[6N] = -J— [ d^x6N[n'^'^ + 167rG'(H« +H(^))] (39) 
IottG J 

with the Hamiltonian densities given by 

n^^^ = -Ak^5]5Ki - ^6i6E] + -^d^d^dE;; (40) 

and 

^1'^ = -37? - ^v=o, H«=p^/^(^K.(^)5v^ + nv^)^j • (41) 

For the diffeomorphism constraint (jlj), we have the second order term 

D^\bX'^\ = ^J^d'x6N'^ [pdM5K',)-p{d,5Kl) - k5t{da5Et) + MG{ixdM)\ ■ 

(42) 

In what follows, we perform canonical gauge transformations on the basic phase variables 
and demonstrate that this, too, results in the gauge invariant combinations (1351) . 
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3.1.1 Gauge transformations of basic variables 

We start by computing gauge transformations of the basic phase space variables: K^, Ef, 
{p and TT. Using (l38l) with the constraints (|39l) and (H2i) we obtain 



^2' 



%„,^]57r = 7rAe-pV'(^)eo (43) 

for the basic gravity and matter perturbations. Note that the background lapse function 
has again been set to iV = y/p for gauge transformations of a metric in conformal time. It is 
also easy to see that when acting upon the background quantities A;, p, and vf these gauge 
transformations will generate only second order contributions. Hence in the equations of 
motion the background phase space variables can be treated as gauge-invariant up to the 
desired order. 

The perturbative treatment introduces a subtlety in the interpretation of transfor- 
mations: In the unperturbed case, the gauge transformation of a phase space function 
X{K, E, if, tt) generated by the diffeomorphism constraint acts as a Lie derivative 

{X,D[C]} = C^^X (44) 

along the vector field At the same time, the diffeomorphism transformation of the 
perturbations is given by the ^-part of fHSl) . i.e. {6X, D[^°-]} = 6[o^^]6X, whereas the barred 
quantities remain intact at the linear level {X,Z)[^°]} = 0(2). For a scalar field, these 
Poisson brackets can directly be identified as Lie derivatives of background quantities and 
perturbations: 

C^^ = Cda^ = , C^^ip = CdaS^ = 0(2) . 

However, if one computes the Lie derivatives of perturbative terms of a tensorial object, or 
even of a scalar of non-zero density weight, one can notice that they do not coincide with 
the ^-part of (H3l) . For instance, 

while 

5[o,5]7r = , S[o^^]Sn = ttA^ = ndaC 

for C,"" = <9"^. Similar discrepancies occur for the triad and extrinsic curvature. 

Nevertheless, gauge transformations are related to the Lie derivative. As shown in 
App. A of up to a higher perturbative order we have 

{X, Die]} = [C,Xf\ {SX, D^'^IC]} = (45) 
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for the background and perturbed parts of X respectively. (As explained in more detail in 
[2j, the second Poisson bracket takes into account the fact that perturbations 6X do not 
contain zero modes.) These two equations above are consistent with the Lie derivative of 
the full variable X 

{X + 6X, D[C]} = C^X + 6X). 

Individual terms in this expansion, when computed either via (jHj) or via (138!) . do not 
agree in general because may not equal C^X^^^) = C^X: While the ^-gauge 

transformation of an unperturbed variable is equivalent to taking a Lie derivative, the 
perturbation procedure breaks this equivalence. For instance, the Lie derivative of a back- 
ground quantity (being linear) contributes to the gauge-transformation of its perturbation, 
not of the background quantity itself. At the same time, the Lie derivative of a linear 
perturbation (being at least a quadratic quantity and hence neglected here) contributes to 
the back-reaction on the background. Their zero- and first-order parts of the Lie derivative 
do contribute to the diffeomorphism transformation, with combined contribution equal to 
the diffeomorphism transformation of the full variable, although in a rather mixed way. 
When combined to X + SX, conventional transformations are obtained. 

3.1.2 Transformation of the lapse function and shift vector 

In the covariant formulation, the lapse function and shift vector are merely components of 
the space-time metric, and hence subject to coordinate (gauge) transformations in the same 
way as any other metric component. In the Hamiltonian framework, on the other hand, 
lapse and shift act as Lagrange multipliers and are not phase space variables. Therefore, 
unlike e.g. triad components, their gauge transformations cannot be directly obtained as 
Poisson brackets fl38l) (which would always give zero) with the gauge-generating constraints. 
Nevertheless, there exists an indirect procedureJl 

A coordinate change causes a change in the space-time foliation by spatial slices, and 
thus the induced and i^* change according to their gauge transformations. Since the 
slicing is determined by lapse and shift they, too, must change. Lapse and shift not only 
determine the slicing but also, as in (!26|) and (128!) . equations of motion which triad and 
extrinsic curvature have to satisfy as one moves from one slice to the next. Consistency 
of equations of motion for the gauge-transformed canonical variables thus requires certain 
transformations of the lapse and shift: They have to change such that they generate the 
correct equations of motion for the transformed and K^^, on which a canonical gauge 
transformation has been applied. In this way, gauge transformations for N and N°- result 
unambiguously, even though these are not phase space variables. 

Hamilton's equations such as fl26l) have a time derivative of a phase space variable 
on their left hand side, while the right hand side depends on phase space variables and 
the Lagrange multipliers. Thus performing a gauge transformation on both sides of the 
equations one can obtain the transformations of 6N and 6N"'. The non-trivial part of this 

^Alternatively, Poisson brackets can be defined on an extended phase space which also includes the 
Lagrange multipliers ^41j . 
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recipe is the gauge transformation of the left hand side, as a time derivative of a phase space 
variable is not itself a phase space variable, and hence its Poisson bracket with constraints 
is not defined. Furthermore, a gauge transformation does not, in general, commute with 
taking time derivatives, that is the gauge transformation of a time derivative is not merely 
given by the time derivative of the gauge transformation. 

Nevertheless, the gauge transformation of a time derivative can be computed with the 
help of 

Lemma 1 For an arbitrary linear phase space function 6X , the commutator between its 
gauge transformation and its time derivative is given (up to second order terms) by a single 
gauge transformation 

%„^](5X)-(%„^]5X)-=5[o,5„]5X (46) 

Proof: Using the definition of the gauge transformation fl38|) and time derivative (l26l) via 
the Poisson bracket, by virtue of the Jacobi identity we obtain 



{{6X,H[N] + D[N-]},H[mo]+D[d'^^]} 
- {{5X,Hmo] + D[d''^]},H[N] + D[N^]}} 
{5X,{H[N] + D[N%Hmo]+D[d^^]}} . 



(47) 



The inner Poisson bracket can be computed using the constraint algebra. In perturbative 
form, we have [2] 



{H[N]+D[N%H[N^o]+D[d''^]} = D 



P 



H[d''idaN] 



Most of these constraint terms are at least of second order, and constraints whose Lagrange 
multiplier is quadratic will not affect the (leading) linear part of the gauge transformation 
of 5X. Therefore the only relevant contribution comes from the first part of the first term. 



D 



N 



9"(iVeo) 



which is equivalent to a single diffeomorphism transformation with the 



(linear part of the) shift vector given by 

-9'^(iVeo)=9% + 0(2). 
p 

The commutator then reads 



which is (US]). □ 
The last equation implies that diffeomorphism invariant canonical variables do have com- 
muting time derivative and gauge transformation. Moreover, the leading diffeomorphism 
term originates from the Poisson bracket of the two Hamiltonian constraints on the left 



16 



hand side of (jUj). Therefore, taking a time derivative commutes (up to quadratic terms) 
with a diffeomorphism transformation. This can also be seen from the absence of ^ on the 
right hand side of ( l46l) . In other words, gauge transformations which do not involve the 
Hamiltonian constraint (i.e. such that = 0) commute with taking time derivatives. 

For later convenience we write out gauge transformed time derivatives for a number of 
phase space variables: 

kod'P) = (%o,c]V')' 
kodE) = (%o,c]^)' + eo 

kod^^) = (%„?]'5vr)-+7rAeo, (49) 

where the gauge transformations of the triad components, following from fHSj) by compar- 
ison with (Q, are given by 

%,5]V^ = -^eo , ko,^]E = C (50) 

We now have all the ingredients to obtain the transformations of lapse and shift perturba- 
tions which (for the scalar mode and in conformal time) are expressed as in (11 01) . Writing 
the perturbation of extrinsic curvature using the equation of motion 

for the perturbed triad and the expression = a/a = p/2p for the Hubble parameter, we 
obtain 



SKi = -5i. 



^ + J^{^ + (j)) +dad' J^E-{B-E) . (51) 



Performing a gauge transformation of the left hand side according to fll5]) and comparing 
to the gauge transformation of the right hand side using (H9|l and (!50|) yields the desired 
transformation of the Lagrange multipliers. Specifically, the diagonal part provides the 
transformed lapse perturbation 

ko,i]<P = ^o + ^^o, (52) 
whereas the off-diagonal part implies 

5K„,e](i?-i^) = -^o. (53) 

Thus, the transformations of lapse and shift are indeed determined by gauge transforma- 
tions of the phase space variables through the dynamical equations of motion. 

Gauge invariant combinations are then obtained from the gauge transformations, which 
reproduces the Bardeen variables fl35l) . Note that, as follows from fl53l) . the quantity 
B — E is diffeomorphism invariant. According to Lemma [1], the gauge transformation of 
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its time derivative is then given simply by the time derivative of its gauge transformation. 
Consequently, the last term in the ^-equation (l35l) gauge-transforms as 



[{B - E)j = -E)\= -L 

which along with the second term of $ compensates the gauge-dependence of 0. 

We have determined gauge transformations of lapse and shift by making sure that the 
form of equations of motion for E!^ and K^^ is invariant. Here, we used the basic fact 
that changing the lapse function and shift vector leads to a different space-time metric 
decomposition and hence a different form of the triad and extrinsic curvature as well as 
their evolution. A different slicing of space-time also affects the matter variables, e.g. 
the definition of the field momentum. Gauge transformations of the matter variables will 
induce transformations of and A^" through Hamilton's equations for and vr as they did 
for gravitational phase space variables. We must therefore ensure that the transformations 
of the lapse and shift generated in this way be consistent with those obtained in ( l52l) and 
( 153|) . Taking the gauge transformation of the left hand side of the equation of motion 

N ASEt6:)\ , 6N 



p3/2 



= — + -(3V-AE + 0) (54) 
p p 

according to ( H9|) and ( l43l) . and comparing it with the transformation of the right hand 
side, yields the transformation of the lapse perturbation, 5[^o,C]0 = + which agrees 
with Eq. (152|) . Repeating the procedure for the momentum equation 



= pA6if -p^V^^6if -V ^{3^- AE + <!))) + fcAB, (55) 
results in S[^q^^]B = ^. The latter along with (H9|) and (!50|) . implying 

%o,C]^= (%o,?]^)'+eo = e + eo, 

reproduces the correct gauge transformation (!53|) . Thus a fixed transformation of lapse 
and shift provides the correct gauge transformation of both gravity and matter phase 
space variables. This is a further consistency property ensured by the first class nature of 
constraints. 

For matter fields, the gauge invariant density and scalar field perturbations are 

5pGi ^^p^ {B-E)- 6^"^' = 5^ + ip{B-E). (56) 

It is convenient for cosmological applications to introduce the gauge-invariant quanti- 
ties: 



n = ^ + je [^^j =tlj + j^ \^-^] (57) 
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and 



{^ = ^ + je =^P + J^^ (58) 

ff)-(T)— (f)-(l) 

The following points are worth noting regarding the above three gauge-invariant quantities: 
(i) TZ provides information about the nature of the long wavelength perturbations i.e. when 
the perturbations have left the Hubble radius. More precisely, IZ vanishes in the long 
wavelength limit if the perturbations are adiabatic [121 113] • (ii) C refers to the 3-curvature 
perturbations on uniform density hypersurfaces. As IZ, C, is also conserved in the large scales 
and quantifies the large angular scale temperature anisotropies in the cosmic microwave 
background. In the slow- roll limit, the two gauge- invariant quantities are identical and 
either of them can be used to quantify the primordial perturbations, (iii) Unlike C, and 
7^, the Bardeen potential $ evolves in time from Hubble exit until the re-entry during 
matter /radiation era. More precisely, during inflation, at the super- Hubble scales, 

^-t^X =^ ^-esM\ (60) 

where e^^ is the slow- roll parameter which is much less than unity and \A\ is the value 
of C at the super-Hubble scales. However, at horizon re-entry, $ ~ ^. (iv) IZ (and also 
Q is linearly related to the Mukhanov- Sasaki variable Q which is useful for studying the 
quantization of perturbations, (v) Unlike TZ and C, "^2 is not often used in the study 
of cosmological perturbations, (vi) In the quantum corrected version of the perturbation 
equations, all the above quantities acquire non-trivial quantum corrections which we will 
discuss in the rest of the paper. This also affects the conservation of power at large scales. 

For the classical gauge transformations we thus produce the well-known gauge invariant 
quantities fl35|) . but we have now done so in a way which is entirely canonical. These 
methods therefore generalize directly to the case of equations and gauge transformations 
which are corrected by effects from quantum gravity even in quantum regimes where no 
underlying smooth space-time picture exists. 



3.2 Inclusion of quantum corrections 

With quantum gravity effects, both the equations of motion and gauge transformations are 
governed by the quantum corrected Hamiltonian constraint H'^ in f|T^ . including all the 
counter-terms and the diffeomorphism constraint fH2|) which remains unaffected. Recall 
also that the terms a'p, v'p as well as the counter-term functions are leading order quantum 
effects (not to be confused with perturbative order). From now on we neglect all higher 
order quantum corrections, such as (a')^p^, a'pf etc. Moreover, also in the background 
equations (l30l) . (13T]) and (!32|) the quantum corrected terms must be used for consistency. 

Gauge transformations are generated by the quantum corrected constraints with the 
lapse function and shift vector parameterized by the infinitesimal vector field = {C,o, 8°"^,) 
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as 

SN = N^o, SN" = d^C (61) 

The gauge transformation of the triad perturbation taking into account the counter-terms 
can be computed using the Poisson bracket 

^ {6E^,HQ^'\m,] + Z}(2)[5a^]} = 2akpUl + m+p{6tA^ - d^d^O (62) 

from which the transformations for ip and E follow 

%,^]^ = -^eo(l + /), ho,&E = ^- (63) 
The gauge transformed extrinsic curvature yields 



%„,5](aif:) =9U(^e + «%)- 



+ g)- AnG^oaipVi^) - ^(1 + /s)) 



61 



(64) 

We compare this equation with the gauge transformation of (127|) . eliminating the potential 
term using the background Raychaudhuri equation fl3ip . By virtue of the anomaly- freedom 
condition (1221) along with [2j 

/2 = 2/i, (65) 

the off-diagonal part results in 

Si^,,^]iB-E) = -a%, (66) 

whereas the diagonal part yields 

= to + mo (l + V'P + ^) • (67) 



Remarkably, the anomaly cancellation condition (]20l) implies that the last two terms inside 
the parenthesis mutually cancel, hence quantum corrected transformation of the lapse per- 
turbation is equivalent to the classical one fl52l) . Finally the matter perturbation transforms 
as 

%„,5]5^ = ^(l + /i)eo. (68) 

The four metric perturbations can be combined into two gauge invariant quantum corrected 
potentials 

^ = ^-^il + f)?^ (69) 

. , B-E\ ^^B-E 
$ = 0+ +^ 
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Similarly, the gauge invariant matter variables are 



5cp''' = 5cp + ^{l + h)^—i^ (70) 



B-E 



Note that when omitting the quantum corrections in Eqs. (169!) and (!70|) one recovers the 
classical results (135|1 . From the corrected gauge invariant expressions, one can directly see 
that the combination 

^ , 1 + f ^ ^ 1 + f ^ , 

which does not refer to the non-trace perturbations E and B, is gauge invariant. Also 
the explicit a-dependence drops out, showing that this particular perturbation is quantum 
corrected only because we were required to include counter-terms / and /i in addition to 
the primary correction function a. Similarly the other curvature perturbation is 

which does not refer to a, either. (A generalization of the perturbation ( to quantum 
corrected equations requires a general derivation of Sp*^^, which will be done elsewhere.) 

4 Gauge invariant equations of motion 

We can now formulate the perturbed equations of motion purely in terms of the gauge 
invariant variables derived in the previous section. The following auxiliary relations will 
be useful: 

6E^ = -2p^6^ 

-2J^p{l + f)^^6t+p{5^A - d^d,)E (73) 
a5K = -5:(v& + ^(^ + $(l + /)) 

-^a^^ (^(1 + /) - + {^^ -{B-E) 

where the first line of each equation contains only gauge invariant terms. 

4.1 Diffeomorphism constraint equation 

Varying the diffeomorphism constraint (j4]) with respect to the shift perturbation yields the 
diffeomorphism constraint equation (the space-time component of Einstein's equation): 

= STrGa^l^^ = p {d,{a6K'Jf) - dk{a6K^)) - akd^SE^S'^ + 8nGa7TdJ<f . (74) 
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Using the gauge invariant variables defined in (169|) and (!70|) . it can be rewritten as 



^ + + /)$ - A7rG-ip6ip^^ + (gauge terms) = 0, 

V 



(75) 



where the 'gauge terms' are 
B-E 



a 



^(1 + /) + ^2 ( 1 + / + ^) - 47rG-0'(l + /i) 



(76) 



After eliminating the Jif-term using the background Raychaudhuri equation flSTl) the ex- 
pression inside the square brackets of (176!) becomes 



which is proportional to one of the anomaly-freedom conditions fl22l) . Thus all gauge 
dependent terms vanish and the diffeomorphism constraint equation takes the form 



(77) 



Note that classically the right hand side is nothing but the gauge invariant space-time 



component of the perturbed matter stress-energy tensor —ArcGa'^STT 



4.2 Hamiltonian constraint equation 

As seen explicitly for the diffeomorphism constraint, for constraints which are part of a 
closed system, i.e. which result from an anomaly-free quantization, gauge invariance of the 
equations of motion is guaranteed and showing that the gauge dependent terms of a given 
equation do vanish is in general a rather tedious, although straightforward, exercise. We 
leave out such explicit demonstrations in this and the following sections. 

The Hamiltonian constraint equation is obtained by variation with respect to the lapse 
perturbation: 



S{SN) 



WttG 



p 




,3/2 



h 



/2 



2(y 



SEfSi 



2p 



(7J 



where 



can be obtained from fll42p . 



2p3/2 



(79) 
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Multiplying both sides by (^l\fp again allows one to replace the background extrinsic 
curvature with the Hubble rate. Then eliminating the field momentum and its perturbation 
using (|29l) and the auxiliary expressions (|73l) along with anomaly-freedom conditions of [2] 
to reduce the number of counter-term functions, one arrives at the gauge invariant Hamilton 
constraint equation (or perturbed Friedmann equation) 



A(a2$) - 3jr(l + /) ^ + ^$(1 + /) 



47rG'-(l + /3) ^(5v9«^-^'(l + /i)$ 
V L 

+ z7p\/^((^)5(^^i] (80) 



Again, the right hand side is nothing but the time-time component of the perturbed stress- 
energy tensor, which now includes quantum corrections. 



4.3 Hamilton's equations 

As mentioned before, each pair of Hamilton's equations for configuration variables and 
momenta can be combined into a single second order equation. We illustrate the procedure 
starting with the matter field. The time derivative of the momentum perturbation can be 
first computed using the Poisson bracket 



+N 



and then compared to the time derivative of the right hand side of Eq. fl2I?l) . which will 
include second time derivatives of the scalar field. With the help of the background equa- 
tions and anomaly cancellations conditions, the Klein-Gordon equation can be cast in the 
gauge invariant form 



5¥3«i + 2e^50^' ^'^^) ~ ''''^^ ~ /a)^'^^^' + ''py,^^ (82) 

+ 2upV,^ (<^)(1 + /!)$ - ^ [(1 + /i)<i> + 3(1 + gi)^] - 2^<^(/^p)$ = . 

Similarly one arrives at the spatial components of Einstein's equation. Taking the time 
derivative of Eq. (!27|) and noting that 

{a6Ki) = a6ki + dK^di = adK^ + 2M'a5Kl (^^^ , 

one can substitute the time derivative of the extrinsic curvature perturbation using the 
Poisson bracket 

5k\ = {5KlH[N] + D[N']}. 

The resulting expression will contain second order correction functions a*^^-* and v^'^^ related 
to the background ones, a and z/, by the conditions (l23l) for anomaly freedom. 
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Taking the trace of each equation and substituting it back into the left hand side, one 
obtains 



a'5E: 



d{5Ef) Qp 
whose gauge invariant parts read 



^ {SiS: - 25151) 



GI 



a' 



-^51 



WW) 



u'5E\ 



Qp 



- mi) 



(83) 



GI 



v' 



3 " 



B4) 



The combined second order equation naturally decouples into two independent equations 
by taking diagonal and off-diagonal terms. After a tedious but rather straightforward 
computation, taking into account the background equations of motion together with the 
anomaly-freedom conditions, the former equation takes the form 



2VI/ ( 1 - ^ ) + $(1 + /) 



^ + 2^'' ( 1 + /'p - ^ 

a 



$(1 + /) 

47rG^ [050^1 - vW^^{ip)5i{P'-\ . (85) 



In the absence of anisotropic stress in the matter sector, which is the case for the scalar 
field, the gauge invariant part of the off-diagonal equation reads: 



a„a^(a2($-^(i + /i))) = o. 



(86) 



which implies $ = \1'(1 -|- /i), replacing the classical relation $ = Here, /i is a counter- 
term correction which has to satisfy 



^7) 



5 Qualitative properties of the scalar perturbations 

In this section, we discuss salient properties of the scalar perturbations including inverse 
triad corrections from loop quantum gravity. For now, we do not compute the power- 
spectrum since this would involve reducing the perturbation equations ( 1771) . ( IHOl) . ( l85i) . 
( 182|) into a single differential equation in terms of the Mukhanov- Sasaki variable. Instead, 
we will focus on aspects of the classical matter perturbations in an effective quantum space- 
time. In particular, we show that (i) the speed of scalar perturbations is less than unity 
and can in fact be much smaller, and (ii) the scalar perturbations are not purely adiabatic 
and have a small entropic contribution arising from the quantum correction. 
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5.1 Speed of perturbations 

Using the relation $ = (1 + the perturbation equations (ffTl) . (ISOl) . (j85l) lead to 
+ /i)A^ - 3jr(l + /) + ^(1 + /)(1 + h)-^ 



GI 



^ + e^(l + /)(1 + /i)^ = 



GI 



2 1- 



a 



+ 



* + + + 



+ 2^^ 1 + /'p - — + 



a 



(89) 
(90) 



Dividing (1551) by 1 + /s and subtracting the result from Eq. (M?]) leads to 



1 + /3 



2(1- — ]+(l + /)(l + /^) + 3^ + ^ 



a 



1 + /i 



+ ^ + [l + fp- 



a'p 
a 



I + /3. 

1 + f){l + h)+J^\l + h) 



(91) 



X 



(1+/)^ i+«w« ,, , ,a , ^ {i+h){i+h) 



^ = -87rGap(l + /a)^^^^^ 



GI 



Replacing V^cp and 5^? using the Klein-Gordon equation ( !32l) and the diffeomorphism con- 
straint fl89l) . we obtain a second-order partial differential equation for 



2^ <^ 1 + - 

I 2 1 + /3 



3/ -/s / 



l + h (p 



2^ Wi + /) + jr 



a 

[1 - v'p/{?,9)] 



(92) 



+ (1 + /) ^ + jr^ <^ 3 



, (1 + /)^ 
1 + /3 



a'p 



+2 1 + - — 1 + /) - -. 37377— ttI 



l + /i)-4 (1 + /) 



where the second term shows the speed of perturbations 



2 -2 

c: = a 



l + h 



(93) 



Eqs. ( !92|) and (!93|) constitute one of the main results of this paper, which has the following 
implications: 

(i) In the perturbative regime analyzed here, we have a > 1 such that there is a danger 
of the speed of sound becoming super-luminal. However, as in the context of gravita- 
tional waves in [20], one must compare the propagation speed not with the classical 



25 



speed of light (which is one) but with the physical speed of light of electromagnetic 
fields in the same effective quantum space-time. Since the Maxwell Hamiltonian is 
subject to quantum gravity corrections, too, [3, |44j the physical speed of light can 
differ from the classical one. In fact, for an anomaly- free coupling of the Maxwell field 
to gravity it must be larger than one by a factor which equals as shown in [20] • The 
speed of perturbations derived here is super-luminal compared to the physical speed 
of light only if the remaining factor (1 + + /a) in fl93|) is larger than one. Using 
Eqs. f l87|) and fl79|) it is easy to show that, in the perturbative regime considered here 
(where a' is negative), we have (l-|-/i)/(l-|-/3) ~ — /s < 1. Thus, in the pertur- 
bative regime, the speed of scalar perturbations is indeed less than unity. Again, one 
can see the importance of consistency conditions for the counter-terms. (Although 
such a scenario arises in the case of non-canonical scalar field inflationary models [15] , 
classically it is not possible for canonical scalar field inflation. Holonomy corrections, 
which have been used in |Tlj without ensuring consistency and anomaly-freedom, 
have been claimed to lead to a speed of sound much larger than unity and even 
divergent in some phases. This may indicate the inconsistency of the perturbation 
equations used there.) 



(ii) From the corrected diffeomorphism constraint equation (1891) one learns that in the 
absence of matter fields, the metric perturbation decays more slowly compared to the 
classical case. Assuming that a is a slowly varying function, the metric perturbations 
decay as 

According to (1571) . f -\-h = 2a'p/a < (and fh is subdominant compared to f + h in 
the regime considered here). Hence, the decay should happen more slowly, implying 
that the inverse triad corrections enhance metric perturbations. 

(iii) In the long wavelength limit, the second term in Eq. fl52]) can be neglected and the 
perturbations can be treated to be independent of the wave-number |k|. Assuming 
that a is a slowly varying function, the Bardeen potential is given by 

3 1 

\1> oc e^^a3-2nc. 2 (95) 



where the constant of proportionality is determined by the choice of quantum state 
defined at the initial epoch of inflation [32]. Comparing this expression with the 
corresponding classical equation (160!) suggests that the primordial perturbations have 
a different behavior compared to their classical counterpart. For < ria < 3/2, 
the quantum perturbations enhance the primordial perturbations compared to the 
classical one. For > 3/2, the perturbations decay and will lead to tiny primordial 
perturbations. This suggests that quantifying the primordial perturbations can, in 
principle, constrain the value of tIq, > 3/2. We will discuss the implications of this in 
a future publication. 
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(iv) Note that in arriving at the above perturbation equation, we have not used the per- 
turbed Klein-Gordon equation (!82l) . To obtain the primordial power spectrum, we 
would need to consider the combined evolution of the scalar and the metric pertur- 
bations 5({)^^ and ^ . 



5.2 Isocurvature perturbations 



For canonical scalar fields in the classical theory, TZ is conserved on large scales implying 
that the perturbations are adiabatic. It can be shown [42] that 



class 







(96) 



Thus, on these scales 7^ is conserved. However quantum effects modify this and lead to 
small entropic perturbations. Taking a time derivative of Eq. ( 1571) and using the diffeo- 
morphism constraint (l89i) leads to 



7^ 



1 + / y 



I jr(l + /)(l + /^) + ^-2| [VP 



(97) 



+2(1 + /)(! + h) - * + ^ ((1 + /)(1 + h))- ^ 

+ 1 ;!^ff_ Vf^ + ^(i + /)(i + /^)^^ 



Eqs. 



(1 + /i)ay V ^ " - ^ , J ^^Q.. 

and fl52]) in the long wavelength limit imply 



n 



l + fia 



a 



I + /3 



^2 



(9J 



from which we infer that TZ is not conserved on large scales. Thus, perturbations generated 
during a single scalar field epoch are no longer purely adiabatic — perturbations contain 
a small entropic contribution. Although such an effect has been seen earlier in a class of 
Lorentz violating models [IB], our example here is the first to show that the primordial 
perturbations from inflation are not purely adiabatic and always contain a small entropic 
perturbation. (Note that anomaly-freedom of our equations [2j ensures that there are no 
violations of Lorentz symmetry, although the speciflc form of symmetries can be quantum 
corrected.) 



6 Discussion 

Effective constraints lead to quantum corrections in equations of motion as well as in 
gauge transformations. Applied to general relativity, it is not only the dynamics but also 
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the underlying space-time structure and the notion of covariance which are affected by 
quantum corrections. In [2] it was shown that the full constraint algebra, i.e. the set of 
structure functions, of canonical quantum gravity changes when quantum corrections of 
a loop quantization are included. While no gauge freedom is destroyed in the anomaly- 
free quantization used, the algebra is not the classical hypersurface deformation algebra as 
originally derived by Dirac for classical general relativity. The quantum corrected equa- 
tions used are generally covariant, but the symmetry type of the underlying covariance is 
quantum corrected. A determination of the full gauge algebra of quantum gravity would 
require going beyond the leading perturbative order, which is not available so far. But the 
results of [2j show that quantum corrections in the algebra must arise. In particular, this 
implies that terms in an effective action of canonical quantum gravity cannot be simply 
of higher curvature form. There must be additional effects such as non-local terms or 
non-commutative manifold structures. 

In [2j as well as in this paper only the inverse triad type of quantum corrections due 
to the effects of [SI E] is considered. There are additional quantum corrections, one due 
to the use of holonomies and a generic one due to quantum back-reaction of fluctuations, 
correlations and higher moments of a quantum state [17]. For these corrections no anomaly- 
free version for perturbative inhomogeneities has been found yet, which indicates that 
there are severe consistency restrictions especially for holonomy corrections. (Quantum 
back-reaction is generic, such that the existence of consistent deformations is guaranteed 
by the work on effective gravity e.g. in [HlllH].) Despite this incomplete status of quantum 
gravity corrections, the different structures of the three types of corrections shows that it 
is not possible to cancel corrections from one type, such as the inverse triad corrections 
used here, by corrections of the other types. Thus, corrections of the form discussed here 
must be present in any cosmological perturbation theory based on loop quantum gravity. 

Given anomaly-freedom of the constraints, it is possible to construct gauge invariant 
variables and recast the equations of motion in an entirely gauge invariant manner. Here, 
this analysis was done for the perturbative constraints of p], incorporating inverse triad 
corrections of loop quantum gravity in a way which, to leading orders, is anomaly-free. 
The final equations (1771). (I5U]). (1521). (|g^. and 
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+2PpV,^ + - ^ [(1 + /!)$ + 3(1 + g,)^\ - 2^(^(/^p)$ = 

collected here from the last two sections are manifestly gauge invariant and reproduce 
the classical perturbed Einstein's equation if one omits the quantum corrections. Besides 
gauge invariance, there is also a consistency issue which arises since, on general grounds, 
there are three unknown scalar functions subject to five equations. Moreover Eqs. fl99|) and 
(11021) can be used to eliminate two of these functions in terms of just one, say $, which 
should satisfy the three remaining equations. 

Another perspective on closure of the equations of motion is given by considering that 
the Klein-Gordon equation (11031) is not independent. In the covariant formalism, it results 
from the energy conservation equation for the matter field: V ^T^'^ = 0, the counterpart 
of the Bianchi identity of the gravitational sector. The latter equation is automatically 
satisfied by construction of the Einstein tensor. For this reason, the Klein-Gordon equation 
can be expressed in terms of the other equations and their derivatives. In the canonical 
formulation such an argument, referring to the Bianchi identity, is not available, especially 
at the effective level, for it is a priori not clear what kind of action might correspond to the 
quantum corrected constraints. Nonetheless, one can use the form of the Bianchi identity 
as guidance to explicitly check the redundancy of the Klein-Gordon equation. In addition, 
the Hamiltonian constraint equation is indeed a constraint, which restricts initial data, 
rather than a dynamical equation as it does not contain second order time derivatives. 
If the constraint is also preserved dynamically it does not break the consistency of the 
equations of motion. 

In the canonical setting, this is realized if the constraints are first class, which is the 
case in the situation at hand to the orders considered. In fact, closure of the constraint 
algebra guarantees both gauge invariance of the equations of motion derived here and their 
consistency. In one of the equations, specifically in the diffeomorphism constraint equation 
f l75|l . we have explicitly shown that all gauge dependent terms mutually cancel. Similar 
straightforward but tedious calculations for the other equations can be performed, but for 
brevity we did not present them here. Given the closure of the constraint algebra, such an 
explicit demonstration of vanishing of the gauge terms becomes unnecessary, although it 
may still serve as an independent consistency check. 

In the light of this, one can arrive at the final gauge invariant equations of motion using 
the following shortcut. After computing the time derivative for the corresponding conjugate 
momentum using the Poisson bracket, it is possible to keep only the gauge invariant parts 
of the variables, dropping all the gauge terms {B, E and {6ip^^ — 6ip) in our case). Although 
setting B = = E in the metric would amount to the longitudinal gauge, this procedure 
is not equivalent to fixing the longitudinal gauge prior to deriving the equations of motion. 
By doing so, one would loose variational equations by off-diagonal metric components 
and thus control on the off-diagonal spatial Einstein equation. Without that equation, the 
relationship between $ and \I' would remain undetermined and one could only 'borrow' the 
relation between $ and from the classical picture. The latter has proven to be incorrect 
at the effective level, as can be seen from Eq. fISBl) . Also effects of counter-terms, required 
for anomaly-freedom, could not be seen in a gauge-fixed analysis, which in general makes 
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such effective equations inconsistent. The same remarks apply to other gauge choices, such 
as uniform gauge. 

Terms which arise from a complete treatment of all gauge properties, but which could 
not be seen in a gauge-fixed analysis, do have physical implications. As an example, we can 
derive an evolution equation for curvature perturbations by subtracting Eq. (IHOl) divided 
by 1 + /a from Eq. fl85|) in such a way that matter perturbations are canceled for a stiff 
fluid presented by the free scalar in the case V{ip) = 0: 

a J 

+ (^^^ + ^M2(l + fp-^)+3^^^-^^l 1$ = 0.(104) 

This equation is similar to fl92|) . but more special because it was derived assuming a stiff 
fluid. For long-wave length modes one can ignore the Laplacian of $ and arrive at an 
ordinary differential equation in time for only if we also use the relation $ = \E'(1 + /i). 
The classical equation would then be solved by a decaying function "^{rj) as well as a 
constant mode because the scale factor for the stiff fluid case satisfies the classical equation 

+ 2J^^ = which makes the coefficient of $ in the evolution equation vanish. This 
conservation of power on large scales |13] can be demonstrated for any perfect fluid by 
eliminating stress-energy components from the gauge-invariant equations. It can also be 
shown to be a direct consequence of the classical conservation law. 

With quantum corrections, however, the coefficient of $ does not cancel exactly and 
the constant mode disappears. First, the background equation is now corrected to 

^ = -2^^(l-iM:^) (105) 

which follows from a combination of fl5Ul) and fl?Il) in the case of a vanishing potential. With 
this, and using our perturbativity assumptions on the correction functions, the coefficient 
of $ in (fTUD is 

^^(2rp+6/-3/3-5^ + |^) (106) 
\ a 3 u J 

which does not have to vanish even if the anomaly cancellation conditions are used. This 
confirms the conclusions of [8j which initially were based on a gauge-fixed treatment in 
longitudinal gauge. However, here the signs of the correction terms are different (for 
instance, /' < while / > 0) such that it depends on the regime whether power is 
enhanced or suppressed. 

The gauge invariant equations of motion can now be used to describe evolution of the 
curvature perturbations, e.g. during cosmological inflation. The small quantum corrections 
accumulated during a sufficiently long inflationary phase may potentially lead to detectable 
imprints on the surface of last scattering and be observed in the Cosmic Microwave Back- 
ground. 
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Appendix 

In this appendix we present the determination of gauge invariant equations of motion, done 
in the main text in the presence of quantum corrections, for the classical canonical theory. 
Resulting equations can be seen to agree with [i2] . 



A Gauge invariant equations of motion 

As we have seen earlier, the background phase space variables as well as the background 
lapse and shift remain unchanged under the infinitesimal gauge transformations (up to the 
second order). Therefore the equations of motion governing them are also gauge invariant. 
In the perturbed context, both the canonical fields and the Lagrange multipliers gauge- 
transform in a non-trivial way, and equations of motion should be formulated for the 
gauge-invariant variables. Here, we systematically derive the background and perturbed 
canonical equations of motion of classical linearized gravity. We also discuss the equivalence 
of the canonical system of equations and Einstein's equation. 



A.l Background equations 



Since the background shift vector is zero, the background diffeomorphism constraint is 
identically satisfied. Therefore background equations are generated only by the background 
Hamiltonian constraint 



H^^\N] = VqN 



' SttG ' 2p3/2 



+ 



7r2 



(107) 



where k and p are the background (diagonal) components of the extrinsic curvature and 
triad respectively, if and vr are the background matter field and its conjugate momentum. 
(The parameter Vq = J d^x is again the coordinate volume of space or of a compact 
region in which the perturbations are introduced.) The background lapse, N = a = y/p, 
corresponds to the conformal time. 

The Hamiltonian (11071) gives rise to the constraint equation 







Vo 



3y^e 

" SttG ^ 2p3/2 



+ p'/V(<^) 
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and two pairs of Hamilton's equations of motion 

k = {^,ffW[iV]} = -^ + 47rGiV(^-^ + Vpn^)j (109) 

p = {p^H^^\N]} = 2N^k (110) 
for the gravitational variables and 



^ = {^,i^(°)[iV]} = ^ (111) 



p3/2 



TT 



{7r,/7W[iV]} = -iVp3/2^^(<^) (112) 



for the matter field. Note that the lapse N = y/p should be fixed consistently after 
computing the Poisson brackets. Then, as follows from Eq. f lllOp . the background extrinsic 
curvature is nothing but the conformal Hubble parameter 

k = f_^-=:J^, (113) 
2p a 

whereas the field momentum is given by 

7f = pLf) . (114) 

Using the two relations above in fllOSp yields the Friedmann equation 

2 



whose right hand side is proportional to the background matter energy density. The 
gravitational equation fll09p should be recognized as the Raychaudhuri equation 

J& = Ji^^ - A-kGi^^ (116) 

written in conformal time. Finally, combining the matter Hamilton's equations into a 
single second order equation results in the Klein-Gordon equation 

^ + 2^^ + pV;^(^) =0. (117) 

The unusual factor of 2 (instead of the standard one, 3) in the second term is again due 
to the choice of conformal time. 



A. 2 Perturbed equations 

Perturbed equations of motion are generated by (the second order part of) both Hamilto- 
nian and diffeomorphism constraints. Using the background and perturbed lapse function, 
the gravitational Hamiltonian constraint can be written as 

H^lm = Y^ld'x [iV7^W + 6Nn^'^] , (118) 
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with the first- and second-order parts of the Hamiltonian density given by (HOj) and 



2k 



-{d,5E'^){dd5El 



respectively. The second order part of the diffeomorphism constraint reads 

1 



In the matter sector, we similarly have the Hamiltonian constraint 



(2) 



matter 



[N] 



N (7^(2) + ^g) + ^(2) ) + + 



;ii9) 



(120) 



;i2i) 



with the densities (HTi) and 



+ ^ 



2p3/2 p3/2 2p 2p3/2 



(2) 



H 



6i6E^ 
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4p 



Yy2 



+p'/'V{^) 



\6i5E^y 5%5E^.5Ei' 



v^2 



along with the diffeomorphism constraint 



^SttejA^I 



4p' 



(122) 



(123) 



In the expressions above, the triad perturbation has the form ([9]), whereas the perturbations 
of lapse and shift are given by ffTOj) . 

Below we formulate the perturbed equations of motion purely in terms of the gauge 
invariant variables. The equations, as before, are of two types: 

(i) Constraint equations, i.e. the Hamiltonian and Diffeomorphism constraints and 

(ii) Dynamical (Hamilton's) equations. 

The latter are those for the matter variables (one second order equation) and for the grav- 
itational variables (two independent second order equations: diagonal and off-diagonal). 
The following auxiliary relation will be useful for deriving gauge-invariant equations: 

5E^= - 2p^5f 

- 2^p{B-E)5''i+p{5'^/\-d''di)E (124) 



5KI 



^ + (^ + <l>) 



5\^{B - E) + dad' J^E-{B-E) 
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where the first hne of each equation contains only gauge invariant terms. 



A. 2.1 Diffeomorphism constraint equation 

Varying the smeared diffeomorphism constraint with respect to the shift perturbation yields 
the diffeomorphism constraint equation (the space-time Einstein equation): 



= SttG 



6D[5N' 



p {dciSK'Jt) - dk{SK^)) - kdd6Et5l + MGiidcd^. (125) 



Using the gauge invariant variables defined in fl35|) and fl56|) . this equation can be rewritten 



as 



^ + - 4:-KGip5ip'^^ + (gauge terms) = 



(126) 



where the 'gauge terms' are 



2{B - E)dc 



+ - 47rG'<^ 



The expression inside the square brackets is nothing but the background Raychaudhuri 
equation f lll6p . Thus all gauge dependent terms vanish and the diffeomorphism constraint 
equation takes the form 



ATtG^pdc^if 



GI 



(127) 



whose right hand side should be recognized as the gauge invariant space-time component 
of the perturbed matter stress-energy tensor —AnGa^STj . 



A. 2. 2 Hamiltonian constraint equation 

The Hamiltonian constraint equation is obtained by variation with respect to the lapse 
perturbation: 



6H 
6{6N) 



1 



WnG 



(— 



0. 



:i28) 



Dividing both sides by ^/p allows one to replace the background extrinsic curvature with 
the Hubble rate. Then eliminating the field momentum and its perturbation in terms of 
the time derivatives of the scalar field (see fl5^ ) and using the auxiliary expressions (11241) . 
one arrives at the gauge invariant Hamilton constraint equation (perturbed Friedmann 
equation) 



47iG 



GI 



(129) 



Again, the right hand side is nothing but the time-time component of the perturbed stress- 
energy tensor. 
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A. 2. 3 Hamilton's equations 



The perturbed dynamical matter equations are computed using the Poisson bracket, giving 
rise to (15^ and (!55l) . Expressing 6it from the first equation and substituting into the second 
one, with the help of the background equations, the Klein-Gordon equation can be cast in 
the gauge invariant form 

6^"^' + 2Jif6^^' - AS^"^' + pV,^^ + 2pV,^ (^)<l> - (^6 + 3^) = . (130) 

Similarly one can arrive at the spatial Einstein equations. The first order equations are 
given by 



ski = {5KIH^^\N]+ D^^^N"]} 
N 



p3/2 



h.2 1.2 Sjik 

-kp5K - -i^5Et5],5l + -5Et5',5l + —dMEt 



1311 



2VF " VP K^Et) 

SE^ = {6E^,H^^\N]+D^^^[N'']} (132) 
= ^ [pSKi6^6'; - p{6Ki6^)6t - dEf] + 25Nk^6t + p i6tdc6N' - d.dN") 

The combined second order equation naturally decouples into two independent equations: 
diagonal and off-diagonal. After a tedious but rather straightforward computation, taking 
into account the background equations of motion, the former equation takes the form 

^ + (2^ + + (^J& + 2^2) $ = AttG (^50^^ - pV^{ip)6^^^) . (133) 

In the absence of anisotropic stress in the matter sector, which is the case for the scalar 
field, the gauge invariant part of the off-diagonal equation reads: 

dad' [$ - ^] = 0, (134) 

which implies $ = \1' for the Bardeen potentials. 



B Anomaly Cancellation Conditions 

In this appendix, we summarize the set of anomaly cancellation conditions containing 
counter-term coefficients. These conditions, in turn, determine the coefficients in terms 
of primary quantum correction functions a, v and a, and they impose restrictions on the 
primary corrections. We note from the expressions of quantum corrected Hamiltonian den- 
sities f lT7|) . ffTSl) that there are three such functions (/, g, h) in the gravitational sector, six 
(/i, /2, 91, 92, 93, 9b) in the kinetic sector and two (/a, 90) in the potential sector of scalar 
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matter. Thus for the system under consideration we have a total of eleven initially unde- 
termined functions contained in the counter-terms. These free functions are constrained 
by anomaly cancellation. 

Invariance of counter-terms under diffeomorphisms pj led to four conditions 

g = -2f , f2 = 2/1 , g, = g,, g, = 2g2 = 2g, (135) 

among these coefficients. These equations trivially lead to the solutions for g, /2, g2 
and gs, leaving seven functions to be determined. Cancellation of anomaly terms from 
the Poisson bracket between Hamiltonian constraints led to two independent conditions 
{Gf = 0, = of [2]) from the gravitational sector. These two conditions along with the 



equation fll35p imply 

h = -/ + 2^ , (136) 
a 

2pf = (137) 

a 

These two equations explicitly solve / and h in terms of the primary correction function 
a. In particular, for the given form of a = 1 + (^p/p)"" + ■ ■ ■, (11371) has the solution 

f=^^ (138) 

2na a 

Thus, there are only five remaining functions that need to be determined. Three of the 
matter anomaly cancellation conditions (Vi = 0, = 0, ^4 = of [2|) lead to 

95 = -fi , 96 = -91 , h = fi- 9i (139) 
whereas another two conditions {B'f = 0, = of [2]) solve /i and gi explicitly as 

/i = /-^, (140) 

91 = ^_El + l(El\p. (141) 
3a u 9 \ u J 

At this stage all counter-terms coefficients have been determined. There is one remaining 
anomaly cancellation condition from the matter sector, P4 = of [2], which implies 

2p/3 + 3/3 - 3/ = (142) 

and thus requires the primary correction functions to satisfy 

a 3 \ a J V 9\v ) 9 \ v ) ' ^ ' 

Independent of the counter-terms and the the requirement (1143p . we also have the relation 
o? = ua to be satisfied by the primary correction functions. Thus, anomaly freedom of 
the constraint algebra severely restricts the allowed form of primary quantum corrections 
functions, but it does permit non-trivial forms of quantum corrections. 
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